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Active  control  of  bending  vibrations  in  a  cantilever  beam 
is  examined  using  a  digital  computer  model  of  beam  and  control¬ 
ler.  The  controller  uses  the  discretized  beam  equation  of 
motion  in  a  linear  control  system,  which  uses  a  Luenberger  ob¬ 
server  to  reconstruct  modal  amplitudes  and  velocities  from  the 
sensor  output.  Feedback  gains  obtained  from  a  steady  state 
optimal  regulator  drive  a  force  actuator.  The  model  is  used 
to  examine  three  areas  of  active  control  of  bending  vibrations. 
First,  impact  on  control  effectiveness  is  investigated  for 
iterative  changes  in  elements  of  the  state  weighting  matrix, 
part  of  the  quadratic  performance  index  minimized  for  the  steady 
state  optimal  regulator.  Second,  the  steady  state  optimal 
regulator  is  replaced  with  classical  control  through  addition 
of  open  loop  zeroes  to  the  system  transfer  function.  Third, 
the  sensor  model  is  changed  to  include  position  and  rate  in¬ 
formation  and  rate  information  only.  State  weighting  matrix 
element  changes  selectively  produce  increased  damping  of  the 
mode  associated  with  the  changed  element.  Breakdown  of  the 
observer  model,  and  instability,  occurs  when  the  change  in  an 
element  exceeds  a  limit  peculiar  to  that  element  and  its  rela¬ 
tive  magnitude.  Control  through  classical  feedback  compensa¬ 
tion  is  at  least  as  effective  as  optimal  control  by  the  steady 
state  optimal  regulator.  Addition  of  rate  information  to  the 
sensor  model  causes  instability  because  of  numerical  inaccura¬ 
cies  in  the  solution  of  the  linear  equation  producing  the 
Luenberger  observer  state  estimation. 
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INTRODUCTION 

Increasing  interest  in  construction  of  large  space 
structures  such  as  solar  energy  collectors  and  large  orbit¬ 
ing  radio  telescopes,  present  unprecedented  problems  in  flexi¬ 
ble  satellite  attitude  control.  Separation  of  control  and 
structural  frequency  bandwidths  to  control  undesirable 
structural  vibrations  may  not  be  possible  for  these  large 
space  vehicles.  Active  controllers  applied  to  structural  vi¬ 
brations  may  provide  an  answer. 

Balas  (Ref  1)  described  a  method  for  applying  active  con¬ 
trol  to  a  simply  supported  beam.  The  equation  of  motion  of 
the  beam  in  bending  is  discretized  by  the  normal  modes  approx¬ 
imation.  This  information  is  incorporated  in  a  linear  control 
system,  where  information  from  the  sensor  is  used  to  estimate 
the  modal  amplitudes  and  velocities  through  a  Luenberger  ob¬ 
server.  The  steady  state  optimal  regulator  produces  control 
proportional  to  the  modal  amplitudes  and  velocities  by  mini¬ 
mizing  a  performance  index.  This  control  is  applied  to  the 
beam  by  a  force  actuator. 

Hungerford  (Ref  4)  implemented  Balas’  method  in  a  digital 
computer  model  of  a  cantilever  beam.  He  showed  that,  at  least 
within  the  limits  of  his  investigation,  the  method  of  control 
appeared  feasible  for  controlling  bending  vibrations  in  a 
cantilever  beam.  Hungerford  looked  at  the  influence  of  con¬ 
trol  weighting,  sensor  and  actuator  location,  and  observer 
error  on  control  effectiveness.  In  this  thesis,  the  influence 
of  state  weighting  changes,  different  sensor  models,  and  the 
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application  of  classical  control  techniques  are  investigated. 

The  steady  state  optimal  regulator  applied  to  a  finite 
number  of  states  in  a  beam  provides  incomplete  control;  that 
is,  an  infinite  number  of  states  remain  uncontrolled  by  the 
regulator  regardless  of  the  number  of  states  controlled. 

For  such  a  system,  there  is  a  need  for  determining  the  effects 
of  changing  the  weighting  in  the  quadratic  performance  index, 
which  is  minimized  to  produce  the  "optimal"  regulator. 
Hungerford  investigated  the  effect  of  changing  the  control 
weighting  matrix  on  the  control  effectiveness.  The  first  pur¬ 
pose  of  this  thesis  was  to  complete  the  study  of  the  effect  of 
weighting  changes  by  determining  the  effect  of  state  weighting 
matrix  element  changes  on  the  effectiveness  of  control. 

During  the  study  of  the  optimal  regulator,  the  root  locus 
indicated  the  effects  of  weighting  changes.  Examination  of 
several  root  locus  diagrams  suggested  the  application  of 
classical  control  techniques.  By  placing  zeroes  at  selected 
points  on  the  root  locus,  it  was  expected  that  effective  con¬ 
trol  might  be  acheived.  Examination  of  simple  feedback  com¬ 
pensation  applied  to  a  vibrating  beam  became  a  second  purpose 
for  this  thesis. 

The  Luenberger  observer  in  Hungerford 's  study  used  input 
from  a  position  sensor.  Foreseeing  application  to  physical 
systems,  some  alternatives  in  selection  of  sensors  is  desirable. 
The  third  purpose  of  this  thesis  was  to  expand  Hungerford' s 
study  to  include  a  rate  sensor  and  a  position-rate  sensor. 

From  a  starting  point  offered  by  Balas,  Hungerford  began 
an  investigation  of  the  active  control  of  bending  vibrations 
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in  a  cantilever  beam.  The  overall  purpose  of  this  thesis  was 
to  further  that  investigation.  The  specific  areas  examined 
were:  the  effect  of  changing  elements  of  the  state  weighting 

matrix;  the  application  of  simple  feedback  compensation;  and 
the  use  of  rate  and  position-rate  sensors. 
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ANALYTIC  TECHNIQUES 


The  basic  vehicle  for  most  of  the  investigation  was  a 
digital  computer  model  of  a  cantilever  beam,  Luenberger  obser' 
ver,  and  optimal  regulator  developed  by  Hungerford. 

For  the  investigation  of  state  weighting  changes  and 
different  sensor  types,  control  was  applied  to  the  first  two 
normal  modes  of  bending  vibration,  with  a  third  mode  included 
as  a  residual  mode.  This  simplified  the  study  without  loss 
of  generality,  and  assured  compatability  with  the  "Total" 
program,  developed  by  Larimer  (Ref  5)  for  control  systems 


analysis  and  design.  "Total"  was  used  to  generate  root  locus 


external  forces: 
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where  the  stiffness,  EI,  and  the  mass  per  unit  length,  m,  are 
constant  along  the  beam  length.  The  equation  is  discretized 
by  using  a  normal  mode  approach  where: 

n 

y  (x,  t)  =  0i(x)ui(t)  (2) 

The  normal  modes  (x)  are  found  by  solution  of  the  special  por¬ 
tion  of  eq  1,  and  they  represent  the  mode  shapes.  The  modal 
amplitudes  u.  (t)  form  the  basis  of  the  linear  system  model. 

The  linear  system  is  constructed  on  a  state  vector  £ 
where : 
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The  subscripts  c  and  r  refer  to  the  controlled  and  residual 
modes,  respectively.  The  linear  system  representing  the 
cantilever  beam,  Luenberger  observer,  and  steady  state  optimal 
regulator  has  the  form: 
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The  subscripts  on  the  control  and  output  coefficient  matrices 
^b]  and  £c]  refer  to  those  coefficients  associated  with  the 
controlled  and  residual  mode,  respectively.  For  a  complete 
discussion  of  the  assembly  of  this  linear  system  matrix,  the 
reader  is  referred  to  Ref  4. 

For  investigating  system  time  response,  eq  4  can  be  in¬ 
tegrated.  However,  a  root  locus  may  be  constructed  from  the 
open  loop  transfer  function  which  has  the  form; 

0  -  fCc]T  [’tl]  -  tAj]  -1  (Bcj  (5) 

The  matrices  appearing  in  this  formula  may  be  obtained  from 
the  construction  of  the  linear  system  model.  The  solution  of 
eq  5  and  the  generation  of  the  associated  root  locus  diagram 
may  be  accomplished  merely  by  loading  the  appropriate  matrices 
into  the  "Total"  program  and  executing  the  required  options. 

An  example  of  the  use  of  "Total"  for  generation  of  the  open 
loop  transfer  function  and  root  locus  diagram  may  be  found  in 
Appendix  A. 

To  investigate  a  certain  state  weighting  configuration 
or  sensor  type,  the  following  steps  were  taken  to  produce  a 
root  locus  diagram  for  analysis  of  control  effectiveness. 
First,  beam  length,  width,  thickness,  stiffness,  and  mass 
per  unit  length  were  selected.  Sensor  and  actuator  locations, 
type  of  sensor,  observer  pole  offset,  and  state  and  control 
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weightings  were  input  to  the  linear  system  model  with  the  beam 
parameters.  The  coefficient  matrix  in  eq  4  was  produced  along 
with  the  component  matrices  [ Ac ] ,  [Ar],  [Bc] ,  |Br},  [Cc} ,  [Cr]  , 
(g|,  and  (kJ.  At  this  point,  if  time  response  information 
was  desired,  the  linear  system  model  was  numerically  integrated 
The  matrices  required  for  solution  to  eq  5  were  loaded  into 
"Total"  and  a  root  locus  diagram  was  generated.  For  the  in¬ 
vestigation  of  classical  control  application,  the  same  steps 
were  followed.  However,  options  of  the  "Total”  program  allow 
replacement  of  the  feedback  transfer  function  produced  by  the 
steady  state  optimal  regulator  with  any  desired  transfer  func- 


EFFECT  OF  CHANGING  ELEMENTS  OF  STATE  WEIGHTING  MATRIX 


There  have  been  several  suggestions  from  various  sources 
for  picking  the  elements  of  the  diagonal  state  weighting  matrix 
[F]  in  the  performance  index: 

J  =  1/2  J  (  v£[F]v,  +  fTRf)dt  (6) 

For  instance,  Bryson  and  Ho  (Ref.  2)  present  a  scheme  for  se¬ 
lection  of  these  elements  based  on  the  estimated  maximum  values 
of  the  states.  Balas  suggests  energy  weighting.  Neither  of 
these,  when  applied  to  the  system  in  preliminary  studies,  yield¬ 
ed  satisfactory  response.  The  first  investigation  in  this 
thesis,  then,  was  a  systematic  exploration  of  the  effect  of 
changing  each  element  of  the  state  weighting  matrix. 

As  a  starting  point,  a  configuration  was  selected  that 
was  previously  shown  to  be  stable  using  the  weightings  as  pick¬ 
ed  by  Hunger ford.  An  aluminum  beam  3  meters  long,  3  cm  wide 
and  .5  cm  thick  was  chosen.  The  observer  poles  were  set  ten 
units  to  the  left  of  the  average  real  value  of  the  system  poles. 
For  the  purpose  of  this  examination,  sensor  and  actuator  loca¬ 
tions  were  picked  sufficiently  far  from  any  node  to  delete  the 
effects  of  locating  close  to  a  node.  On  a  3  meter  cantilever 
beam,  a  position  sensor  at  1.4  m  from  the  clamped  end  and  a 
force  actuator  at  0.6  m  from  the  clamped  end  provides  the  de¬ 
sired  conditions. 

To  begin  the  study,  elements  of  the  state  weighting 
matrix  were  set  to  unity;  that  is 


-8- 


- — ^  ill  iffc'idriilfttt.-  ^  rdt* 


[F] 
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The  elements  were  separately  varied  from  a  value  of  .01  to 
100,  incremented  by  a  factor  of  10.  Analysis  of  the  results 

of  variation  about  a  base  value  of  1  suggested  that  a  more 

.  2 
suitable  base  value  would  be  of  order  10  .  Root  locus  dia¬ 
grams  for  state  weighting  element  variation  from  1,  Figs  Bl 
to  B5  in  Appendix  B,  show  significant  changes  only  when  ele¬ 
ments  are  changed  to  100.  Unless  otherwise  noted,  the  root 
locus  diagrams  and  time  response  results  for  this  and  all  parts 
of  the  investigation  are  for  unity  gain.  The  gain  is  the  co¬ 
efficient  k  in  the  following  equation: 

u'  =  k  (g)  v^  (8) 

where,  again,  u'  is  the  control  input,  (g]  is  the  feedback 
gain  matrix,  and  v^  is  the  controlled  portion  of  the  state 
vector . 

To  provide  a  more  effective  basis  for  the  examination  of 
effect  of  weighting  charges,  a  new  value  was  arbitrarily  se¬ 
lected  as  a  base  from  which  each  diagonal  element  was  alter¬ 
nately  varied.  Since  significant  changes  in  the  root  locus 

2 

were  previously  obtained  with  weightings  of  order  10  ,  a  value 
of  500  was  selected  as  a  new  base  value.  Each  diagonal  ele¬ 
ment  was  set  to  500,  and  in  alternate  executions  of  the  com¬ 
puter  program,  each  element  was  incremented  from  500  to  1,000 
to  2,000.  Again,  a  root  locus  for  each  case  may  be  found  in 
Appendix  B,  Figs  B6  to  B9. 

For  these  weighting  variations,  significant  effects  of 
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weighting  changes  were  noticed  in  the  root  locus  diagrams. 

The  resulting  movement  of  open  loop  zeroes  in  some  case  pro¬ 
duced  changes  in  the  branch  patterns.  More  importantly  the 
closed  loop  pole  locations  for  a  unity  gain  system  shifted 
significantly  with  respect  to  one  another.  For  variations  of 
the  first  two  state  weightings,  these  shifts  were  only  slight¬ 
ly  noticeable.  Augmenting  these  first  two  state  weighting 
elements  to  values  of  10,000  and  100,000  produced  more  notable 
changes  in  the  root  locus  diagrams. 

At  this  point  in  the  study  it  was  realized  that  the  value 
of  the  investigation  was  limited  by  a  major  problem  with  the 
investigative  technique.  Up  to  this  point  the  third  residual 
mode  and  the  Luenberger  observer  were  omitted  when  the  compo¬ 
nent  matrices  were  input  to  "Total”  program.  Although  "Total" 
is  capable  of  handling  the  10  x  10  system  resulting  from  four 
controlled  mode  states,  four  observer  error  states,  and  two 
residual  mode  states,  the  open  loop  transfer  function  produced 
had  a  17th  order  polynomial  numerator  and  an  18th  order  poly¬ 
nomial  denominator.  The  order  of  this  transfer  function  was 
too  large  to  be  processed  by  the  "Total"  program,  so  the  re¬ 
sidual  mode  and  observer  had  been  omitted.  However,  eight 
pole-zero  cancellations  occur  because  the  forward  loop  trans¬ 
fer  function  zeroes  are  cancelled  with  the  feedback  transfer 
function  poles.  Omitting  these  poles  and  zeroes  produced  a 
transfer  function  small  enough  for  the  root  locus  generation 
by  "Total".  In  this  way  the  investigation  was  continued  in¬ 
cluding  the  residual  mode  and  the  observer. 

Figure  1  shows  the  resulting  analysis  performed  on  the 
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3  meter  beam  with  the  first  two  vibration  modes  controlled 
and  the  third  mode  with  no  control  applied.  Each  element  of 
the  diagonal  weighting  matrix  had  a  value  of  500.  The  open 
loop  poles  appear  on  the  root  locus  as  x's.  The  four  open 
loop  poles  grouped  together  near  the  real  axis  represent  the 
observer  open  loop  poles.  The  three  open  loop  poles  located 
on  the  imaginary  axis  represent  the  beam  model  with  no  nat¬ 
ural  damping,  and  their  placement  on  the  imaginary  axis 
represents  the  free  vibration  frequencies  of  each  mode.  The 
small  circles  on  the  diagram  represent  open  loop  zeroes,  and 
the  closed  loop  system  poles  for  a  gain  of  unity  appear  as 
triangles.  It  is  the  relative  positions  of  these  closed  loop 
poles  that  is  of  interest  in  this  investigation.  Changes  in 
the  locations  of  the  closed  loop  poles  on  the  branches  extend¬ 
ing  from  the  imaginary  axis  represent  changes  in  damping  of 
the  associated  mode.  These  changes  are  of  particular  interest 
to  this  investigation. 

The  system  that  appears  in  Fig  1  serves  as  a  basis  for 
comparison.  State  weighting  matrix  changes  will  produce 
changes  in  closed  loop  pole  locations  from  those  in  this  fig¬ 
ure.  These  shifts  represent  changes  in  the  system  response. 
Changes  in  response  can  be  verified  by  examination  of  time 
response  plots,  which  were  obtained  by  numerical  integration 
of  the  system  as  formulated  in  Hungerford's  program.  Time 
response  plots  of  actuator  force,  displacement  at  sensor 
location  modal  amplitudes  of  the  first  three  modes,  and  the 
percentage  of  the  initial  energy  left  in  the  beam  can  be 
found  in  Appendix  B,  Figs  Bll  to  B25,  for  each  system 
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Figure  1:  Root  locus  diagram  for  base  weighting  configuration, 
discussed  in  this  section. 

As  before,  weightings  in  the  diagonal  state  weighting 
matrix  were  changed  alternately,  with  one  element  varied  as 
the  other  three  stayed  fixed  at  a  value  of  500.  The  amount 
of  augmentation  of  each  element  was  based  on  the  previous 
analysis  omitting  the  observer  and  residual  mode. 

A  root  locus  for  a  diagonal  weighting  matrix  with  the 
first  state  weighting  of  100,000  and  the  remaining  weights  at 
500  is  shown  in  Fig  2a.  The  most  obvious  effect  of  this 
heavy  weighting  on  the  first  mode  position  state  is  the 
shift  to  the  left  of  the  closed  loop  pole  on  the  branch 
starting  at  the  first  mode  open  loop  pole.  The  effect  of 
such  shift  on  system  response  is  higher  damping  on  first 
mode  vibrations,  and  examination  of  the  time  response  plots, 

Fig  B15,  verifies  increased  first  mode  damping.  This  weight¬ 
ing  also  causes  changes  in  the  pattern  of  observer  open  loop 
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pole  positions;  all  four  observer  open  loop  poles  are  located 
off  the  real  axis.  This  shift  of  observer  poles  causes  a 
change  in  the  branch  patterns.  The  branch  from  the  third  mode 
open  loop  pole  is  moved  closer  to  the  imaginary  axis,  result¬ 
ing  in  a  shift  to  the  right  of  the  third  mode  closed  loop  pole 
and  reduced  damping  the  third  mode.  Examination  of  the  time 
response  plots.  Fig  B16,  shows  reduced  damping  in  the  third 
mode  for  heavy  weighting  of  the  first  state. 

Second  state  weighting  was  augmented  to  a  value  of 
10,000,  as  shown  in  Fig  2b.  In  this  case,  augmentation  of 
the  weighting  produced  no  noticeable  changes  in  the  closed 
loop  pole  locations,  and  no  noticeable  changes  in  the  time 
response.  The  observer  open  loop  poles,  however,  have  been 
shifted  farther  away  from  each  other.  Augmenting  this  weight¬ 
ing  to  a  value  of  100,000  produces  further  separation  of  these 
observer  poles,  which  causes  a  shift  in  the  third  mode  branch 
into  the  right  half  plane  and  instability  results. 

As  with  increased  first  state  weighting,  increasing  the 
third  state  weighting  produces  a  shift  to  the  left  of  the  first 
mode  closed  loop  pole,  as  in  Fig  3a.  Time  response  plots  for 
the  first  two  modes  are  nearly  identical  for  a  first  state 
weighting  of  100,000  and  the  third  state  weighting  of  2,000. 
This  indicates  that  the  same  behavior  can  be  produced  with 
heavy  weightings  on  position  state  or  relatively  light  weight¬ 
ing  on  the  corresponding  velocity  state.  The  behavior  of  the 
observer  poles  is  not  similar  for  the  two  cases;  for  the  third 
state  weighting  increase,  the  observer  open  loop  pole  locations 
are  not  significantly  shifted.  Consequently,  the  third  mode 
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branch  of  the  root  locus  in  Fig  3a  has  not  been  significantly 
shifted  from  its  position  when  all  weightings  have  a  value  of 
500.  Third  mode  time  response  plots,  Figs  B20,  B21  and  B22 
verify  that  response  of  this  mode  does  not  significantly 
change  when  the  third  state  weighting  is  increased  to  2,000. 

When  the  fourth  state  weighting  is  increased  from  a  value 
of  500,  the  observer  open  loop  poles  shift  off  the  real  axis. 
As  in  previous  cases,  this  causes  a  critical  shift  in  the 
third  mode  branch  pattern,  and  for  unity  gain  an  unstable 
system  results.  To  investigate  the  effect  of  changing  this 
fourth  state  weighting,  the  value  was  decreased  from  the  base 
value  of  500  rather  than  increased.  Figure  3b  shows  a  root 
locus  diagram  for  the  fourth  state  weighting  of  100.  A  shift 
to  the  right  of  the  second  mode  closed  loop  pole  resulted, 
indicating  the  expected  decreased  damping  of  the  second  vi¬ 
bration  mode.  The  time  response  plots  of  this  second  mode. 

Fig  B24,  shows  decreased  damping.  The  shift  of  this  closed 
loop  pole  had  sufficient  numerical  impact  to  cause  a  decrease 
in  the  average  of  the  real  part  of  the  closed  loop  system 
poles  for  unity  gain.  Since  the  observer  open  loop  poles  are 
picked  to  be  ten  units  to  the  left  of  the  average  real  part, 
the  open  loop  observer  poles  are  shifted  to  the  right.  The 
resulting  shift  in  the  branch  patterns  produced  higher  damp¬ 
ing  of  the  third  mode  for  unity  gain,  as  exhibited  in  the 
third  mode  time  response  plot.  Fig  B25. 

The  results  of  this  part  of  the  study  can  be  summarized 
by  examining  the  effects  of  weighting  changes  on  the  closed 
loop  vibration  mode  poles  and  the  open  loop  observer  poles. 
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When  an  element  of  the  state  weighting  matrix  is  increased 
with  respect  to  the  other  elements,  the  associated  closed 
loop  pole  reflects  this  increase  in  a  shift  to  the  left, 
indicating  increased  damping.  For  instance,  relative  in¬ 
creases  in  the  first  and  third  state  weightings  move  the  first 
mode  closed  loop  pole  to  the  left.  Larger  relative  increases 
of  the  position  state  weightings  are  required  to  produce  the 
same  closed  loop  pole  shift  as  a  smaller  relative  velocity 
state  weighting  increase.  There  is  a  numerical  limit  to  the 
tolerable  amount  of  relative  increase  peculiar  to  each  state 
weighting  element  and  the  relative  magnitudes  of  the  weighting 
configuration.  When  this  limit  is  approached  the  mechanism 
for  choosing  open  loop  observer  poles  produces  an  increasingly 
scattered  pattern  of  pole  placement.  This  scattered  pattern 
may  cause  an  outward  shift  of  higher  mode  closed  loop  poles, 
and  instability  may  result  for  increasingly  scattered  patterns. 
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CLASSICAL  FEEDBACK  COMPENSATION  APPLIED  TO  VIBRATING  BEAM 


The  steady  state  optimal  regulator  provides  a  method  for 
controlling  bending  vibration  in  a  beam  optimally  through  mini¬ 
mization  of  a  performance  index.  As  evidenced  in  the  previous 
section,  the  choice  of  weightings  in  this  performance  index 
predetermines  the  limits  of  optimization;  in  some  cases,  the 
optimum  was  still  an  unstable  system.  Although  intelligent 
selection  of  weightings  can  produce  effective  control,  an  al¬ 
ternative  control  law  may  also  effectively  control  bending  vi¬ 
brations  in  a  cantilever  beam. 

In  applying  the  "Total"  root  locus  generator  during  in¬ 
vestigation  of  state  weighting  changes,  it  became  apparent 
that  selection  of  location  of  open  loop  poles  and  zeroes  was 
the  net  effect  of  the  application  of  the  optimal  regulator 
control  law.  The  purpose  of  the  investigation  of  classical 
control  was  to  determine  if  satisfactory  control  of  cantilever 
beam  bending  vibrations  might  be  obtained  by  using  classical 
feedback  techniques  as  an  alternative  to  the  modern  optimal 
control  theory. 

Classical  feedback  compensation  as  applied  to  bending  vi¬ 
brations  in  a  cantilever  beam  deals  with  the  relationships  of 
the  locations  of  system  open  loop  poles  and  zeroes.  By  add¬ 
ing  poles  or  zeroes  at  specified  positions  on  a  system  root 
locus,  branch  patterns  can  be  moved.  The  root  locus  of  the 
system  for  a  3  m  beam  as  described  in  the  previous  section  and 
a  steady  state  optimal  regulator  with  all  state  weightings  set 
to  500  is  shown  in  Fig  4a.  Sensor  and  actuator  positions 


remain  at  1.4  m  and  0.6  m  from  the  clamped  end,  respectively. 
The  Luenberger  observer  has  been  omitted  since  the  states  need 
not  be  known  in  order  to  apply  classical  feedback  compensation. 
Classical  feedback  compensation  was  applied  by  first  removing 
the  steady  state  optimal  regulator  poles  and  zeroes  from  the 
transfer  function  of  the  system  in  Fig  4a.  A  feedback  com¬ 
pensation  transfer  function  was  constructed  to  add  zeroes  in 
the  desired  locations  on  the  root  locus.  This  feedback  trans¬ 
fer  function  replaced  the  steady  state  optimal  regulator  in 
the  linear  system  model. 

The  technique  used  in  this  investigation  involved  the 
use  of  state  space  input,  transfer  function  manipulation,  and 
root  locus  and  time  response  options  of  the  "Total"  program. 
First,  the  system  model  for  the  open  loop  beam  was  produced 
by  the  execution  of  the  digital  computer  linear  system  model 
developed  by  Hungerford.  By  using  the  state-space  input  op¬ 
tions  of  "Total",  the  [A],  [b] ,  and  [c}  matrices  were  loaded 
for  the  uncontrolled  beam.  These  matrices  constitute  a 
6x6  system  of  one  position  state  and  one  velocity  state  for 
each  of  three  modes.  The  forward-loop  transfer  function  was 
generated  from  this  information,  and  the  desired  feedback- 
loop  transfer  function  was  loaded  using  the  transfer  function 
manipulation  options.  The  open-loop  transfer  function  was 
generated  and  used  to  construct  a  root  locus.  Time  response 
to  step  input  plots  generated  by  time  response  options  may  be 
found  in  Appendix  B,  Figs  B26,  B27,  B28,  and  B29.  A  compari¬ 
son  criterion  was  applied  to  these  time  response  plots.  Time 
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Figure  5:  Root  locus  diagrams  for  feedback  conpensation . 


to  settle  to  within  5%  of  the  final  sensor  displacement  for 
a  step  input  was  selected  as  the  comparison  criterion  for  this 
investigation. 

After  examining  Fig  4a,  it  was  decided  that  favorable  re¬ 
sponse  characteristics  might  be  obtained  by  pulling  the  branches 
toward  the  left  of  the  root  locus  after  removing  the  optimal 
regulator  poles  and  zeroes.  Classical  feedback  compensation 
rules  suggest  that  root  locus  branches  may  be  moved  to  the 
left  by  the  addition  of  zeroes  to  the  system  (Ref  3) .  The 
simplest  application  of  the  rule  is  the  addition  of  one  zero. 

On  a  root  locus  diagram,  the  transfer  function  s  +  a  places  a 
zero  on  the  real  axis  at  -a.  Root  locus  diagrams  were  produced 
and  studied  for  location  of  this  zero  at  the  origin  and  2,  4, 

6,  8,  10,  20,  30,  40,  and  50  units  to  the  left  of  the  origin 
on  the  real  axis.  Figure  4b  shows  the  root  locus  for  a  =  0, 
and  Fig  5a  is  a  root  locus  for  a  =  50.  As  expected,  the 
branch  patterns  in  both  cases  are  moved  to  the  left.  However, 
small  values  of  "a"  move  the  lower  frequency  branches  more, 
and  higher  values  of  "a"  have  a  greater  effect  on  the  higher 
frequencies.  Figure  4b  represents  the  former  case,  and  Fig 
5b  represents  the  latter.  Thus,  when  "a"  is  zero  and  the 
gain  is  35,  the  time  to  settle  to  within  5%  of  final  displace¬ 
ment  is  3.2  sec,  as  shown  in  Fig  B26.  The  same  settling  time  is 
3.0  sec  when  the  steady  state  optimal  regulator  is  applied  with 
all  state  weightings  equal  to  500,  as  shown  in  Fig  B27.  When 
"a"  is  50,  there  is  little  first  mode  damping  and  the  5% 
settling  time  is  nearly  infinite,  as  shown  in  Fig  B28. 

It  was  found  that  a  zero  near  the  origin  affected  low 
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frequencies  and  a  zero  far  to  the  left  of  the  origin  affected 
high  frequencies.  An  attempt  was  made  to  affect  both  by  add¬ 
ing  two  zeroes.  One  zero  was  placed  at  the  origin  and  a  second 
zero  was  added  to  produce  the  transfer  function  s(s  +  a). 

Root  locus  diagrams  were  studied  for  location  of  the  second 
zero  at  20,  40,  60,  and  80.  Figure  5b  shows  the  unity  gain 
closed  loop  poles  when  a  =  40,  which  is  representative  of  the 
considered  cases.  The  5%  settling  time  for  this  compensation 
is  2.2  sec,  as  shown  in  Fig  B29. 


RATE  AND  POSITION-RATE  SENSORS 


For  a  cantilever  beam  in  bending  vibration  with  the  first 
three  modes  considered  and  the  first  two  modes  controlled  and 
observed,  the  resulting  linear  system  model,  as  before,  has 
ten  states.  There  are  two  controlled  mode  position  states, 
two  controlled  mode  velocity  states,  four  observer  error 
states,  and  one  residual  mode  position  state  and  velocity 
state.  The  output  matrix  [c]  determines  the  nature  of  the 
output  y;  for  Hungerford's  study  and  the  previous  portions 
of  this  study  the  $c]  matrix  for  a  position  sensor  had  the 
following  form: 

K  c2  0  0  ;  0  0  0  0  :  c5  o]  (12) 

The  more  general  form  of  [c]  is 

f°l  C2  C3  C4  :  0  0  0  0!  c5  c6]  (13) 

For  a  rate  sensor,  c^,  C2»  and  c^  are  zero,  and  for  a  sensor 
providing  a  combination  of  position  and  rate  information,  c^ 
through  c^  are  non-zero.  An  attempt  was  made  to  examine  the 
use  of  a  rate  sensor  and  a  position-rate  sensor  in  the  linear 
system  model  of  the  vibrating  beam. 

Changing  from  a  position  sensor  required  a  substantial 
modification  of  the  linear  system  model  as  compiled  by  Hungerford. 
The  observer  gain  matrix  fK]  was  found  by  setting  the  eigenvalues 
of  the  matrix  (A-(k}{c)]  equal  to  a  desired  set  of  eigenvalues. 

The  resulting  linear  equation  is  simplified  when  c^,  c^  and 
cc  are  zero.  The  more  general  case  requires  the  solution  of 
more  complicated  system  of  linear  equations: 
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In  equations  (14),  r^,  r2 ,  r3,  and  r4  are  the  desired  observer 
eigenvalues.  Solution  of  (14)  using  a  library  linear  equation 
solving  routine  allowed  a  choice  of  sensor  type. 

For  this  part  of  the  study,  a  beam  one  meter  long  was 
used.  All  other  parameters  were  kept  the  same  as  those  in  the 
first  part  of  this  investigation,  with  all  state  weightings 
equal  to  500.  Sensor  and  actuator  positions  were  varied  along 
the  length  of  the  beam  at  0.2  m  intervals,  and  the  linear  sys¬ 
tem  model  compiled  for  each  configuration.  This  was  done  for 
a  position  sensor,  a  rate  sensor,  and  a  position-rate  sensor. 
The  eigenvalues  of  corresponding  systems  were  compared  for 
stability  information.  Computer  model  outputs  for  three  dif¬ 
ferent  sensors  at  0.6  m  and  actuator  at  0.4  m  given  in  Appendix 
C. 

It  was  found  that  including  rate  information  in  the  sensor 
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Figure  6:  Root  locus  diagrams  for  position  and  rate  sensors 


i 


model  caused  unstable  third  mode  vibrations.  This  instability 
results  directly  from  stiffness  and  ill-conditioning  of  the 
linear  system  of  equations  that  is  solved  for  observer  gains. 

A  comparison  of  root  locus  diagrams  for  a  beam  with  sensor  at 
0.6  m  and  actuator  at  0.4  m  with  position  and  rate  sensors 
appears  in  Fig  6.  Although  the  closed  loop  observer  eigen¬ 
values  are  the  same  for  both  systems,  the  figure  shows  a  shift 
in  open  loop  observer  poles  that  causes  a  shift  to  the  right 
of  the  open  loop  zeroes  when  rate  information  is  sensed. 
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CONCLUSION 


1.  Damping  of  a  controlled  mode  can  be  changed  by  making 
a  change  in  the  corresponding  position  or  velocity  state 
weighting  matrix  element  in  the  quadratic  performance  index 

of  the  steady  state  optimal  regulator.  The  rest  of  the  system 
remains  unaffected  unless  the  value  of  this  element  approaches 
the  numerical  limits  of  the  model.  At  this  limit,  the  obser¬ 
ver  model  begins  to  break  down. 

2.  Classical  feedback  compensation  applied  through  addition 
of  open  loop  zeroes  to  the  system  root  locus  can  provide 
effective  control  of  bending  vibrations  in  a  cantilever  beam. 

3.  The  observer  model  is  adversely  affected  when  rate  in¬ 
formation  is  included  in  the  sensor  model,  whether  this 
model  includes  rate  and  position  sensing  or  rate  sensing 
alone. 


1 
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RECOMMENDATIONS 


1.  The  study  of  state  weighting  element  changes  was  limited 
to  a  beam  with  two  controlled  modes  and  one  residual  mode 

by  limitations  of  the  "Total"  program.  The  study  should  be 
expanded  to  consider  a  higher  number  of  modes,  perhaps  by  ex¬ 
panding  the  capability  of  this  package. 

2.  Based  on  the  relationships  discovered  in  this  investiga¬ 
tion  and  in  Hungerford's  work,  control  effectiveness  of  com¬ 
binations  of  state  and  control  weighting,  observer  pole  off- 
sett,  and  sensor  and  actuator  locations  should  be  investigated. 

3.  Classical  feedback  compensation  should  be  applied  to  very 
long  beams.  Hungerford  found  that  the  optimal  regulator  had 
difficulty  with  long  beams.  The  results  of  this  investigation 
suggest  that  classical  control  may  be  more  effective  because 
of  increased  numerical  simplicity. 

4.  Classical  feedback  compensation  should  be  applied  to  a 
beam  model  that  includes  a  higher  number  of  modes  and  higher 
frequency  vibrations.  The  use  of  more  complex  types  of  com¬ 
pensation  may  also  warrant  investigation.  Other  sensor  and 
actuator  locations  than  the  ones  used  in  this  study  should  be 
investigated  using  classical  compensation. 
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APPENDICES 


Appendix  A 


Typical  Use  of  the  AFIT 


"Total"  Package 


The  following  pages  are  a  listing  of  a  typical  session 
with  the  interaction  computer-aided  design  program  for  con¬ 
trol  system  analysis  called  "Total”.  Options  in  the  package 
are  selected  by  a  numerical  input,  and  the  lists  of  options 
used  can  be  found  at  convenient  locations  in  the  listing. 


ATTACH* TOTAL* IU-AFIT 


PFN  IS 
TOTAL 

PF  CYCLE  NO.  =  001 
COMMAND-  TOTAL 

WELCOME  TO  TOTAL — VERSION  1.4 

TYPE  HELP  FOR  INTRO >  TYPE  9V  FOR  NEW  FEATURES  BULLETIN 
OPTION  >  10 


(10-19)  MATRIX  INPUT  OPTIONS 

*10  LIST  OPTIONS 

*11  AMAT  CONTINUOUS  SYSTEM  MATRIX 

*  12  BMAT  CONTINUOUS  INPUT  DISTRIBUTION  MATRIX 

*  13  CMAT  OUTPUT  CONTRIBUTION  MATRIX 

*  14  DMAT  DIRECT  TRANSMISSION  MATRIX 

*  15  '  KMAT  ST.  VAR.  FEEDBACK  MATRIX 

*  16  FMAT  DISCRETE  SYSTEM  MATRIX 

*  17  GMAT  DISCRETE  INPUT  DISTRIBUTION  MATRIX 

*  18  HELP  USER  SET  UP  STATE-SPACE  MODEL  OF  SYSTEM 

*  19  EXPLAIN  USE  OF  ABOVE  MATRICIES 

OPTION  >  18 

IS  THE  SYSTEM  (1)  CONTINUOUS  OR  (2)  DISCRETE?  >  1 


THE  EQUATIONS  YOU  ARE  ABOUT  TO  INPUT  HAVE  THE  FORM: 

XDOT(T)  =  C  AMAT  IX(T>  +  C  BMAT  TIJ(T) 

Y ( T )  =  C  CMAT  DX ( T )  +  C  DMAT  DU(T> 
WHERE  U ( T)  =  GAIN* (  R ( T )  -  C  KMAT  ]X(T)  ) 

AND  X  IS  A  VECTOR  OF  N  STATE  VARIABLES 
U  IS  A  VECTOR  OF  M  INPUTS 

Y  IS  A  VECTOR  OF  L  OUTPUTS 

ENTER  NO.  OF  STATES r I NPUTS r OUTPUTS  >  10ii*l 

ENTER  AMAT  WITH  10  ROWS  AND  10  COLUMNS. 

ENTER  10  ELEMENTS  PER  ROWS 
ROW  1  >0010000000 

ROW  2  >0001000000 

ROW  3  >  -12,002  0  00000000 

ROW  4  >  0  -471.433  00000000 

ROW  5  >  0  0  0  0  -44.73  106.9  1  0  30.16  0 

ROW  6  >0000  10.6  -25.31  0  1  -7.14  0 

ROW  7  >  0  0  0  0  -167.3  371,0  0  0  104.7  0 

ROW  8  >  0  0  0  0  501.7  -1670  0  0  -338.3  0 


ROW  9  >0000000001 


1 


3 


4 


n 


COL 

ROU 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 


0. 

0. 

-12.00 

0. 

0. 

0. 

0. 

0. 

0. 

0. 


0. 

0. 

0. 

-471.4 

0. 

0. 

0. 

0. 

0. 

0. 


1.000 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 


0. 

1.000 

0, 

0. 

0. 

0. 

0. 

0. 

0. 

0. 


0. 

0. 

0. 

0. 

-44.73 

10.60 

-167.3 

501.7 

0. 

0. 


t 


COL  >  6 

ROW 

1  0. 

2  0. 

3  0. 

4  0. 

5  106.9 

6  -25.31 

7  371.0 

8  -1670. 

9  0. 

10  0. 


8 


0. 

0. 

0. 

0. 

1.000 

0. 

0. 

0. 

0. 

0. 


0. 

0. 

0. 

0. 

0. 

1.000 

0. 

0. 

0. 

0. 


0. 

0. 

0. 

0. 

30.16 

-7.140 

104.7 

-338.3 

0. 

-3697. 


10 


0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

1.000 

0. 


ENTER  BMAT  WITH  10  ROUS  AND  1  COLUMNS. 

ENTER  10  ELEMENTS  PER  COLUMN  1 

COLUMN  1  >  0  0  -.1174  .5535  00000  -1.112 


COL  >  1 

ROU 

1  0. 

2  0. 


r 


3 

4 

5 

6 

7 

8 
9 

10 


-.1174 

.5535 


0. 

0. 

0. 

0. 

0. 


ENTER  CMAT  UITH  1  ROUS  AND  10  COLUMNS. 
ENTER  10  ELEMENTS  PER  ROUJ 

ROW  1  >  -.5542  1.324  000000  -.3737  0 


COL  >  1  2 

ROU 

1  -.*,**, an  i  ,tm 


3 


o . 


4 


5 


6 


7 


8 


9 


10 


COL  > 

ROW 

1  0.  0.  0.  -.3737  0. 

r' 


IS  THERE  A  DIRECT- TRANSMISSION  (D  MATRIX)— YES  OR  NO?  >  NO. 
DMAT  SET  TO  1  BY  1  ZERO  MATRIX  (OPTION  79) 

IS  THERE  A  STATE-VARIABLE  FEEDBACK  MATRIX— YES  OR  NO?  >  YES 
ENTER  KMAT  WITH  1  ROUS  AND  10  COLUMNS. 

ENTER  10  ELEMENTS  PER  ROW.’ 

ROW  1  >  10.88  111.1  -45.1  21.8  10.88  111.1  -45.1  21.8  0  0 


COL  >  1  2  3  4  5 

ROW 

1  10.88  111.1  -45.10  21.80  10.88 


COL  >  6  7  8  9  10 

ROW 

1  111.1  -45.10  221.80  0.  0. 


THE  STATE-SPACE  REPRESENTATION  IS  COMPLETE. 
OPTION  >  20  : 


(20-29)  BLOCK  DIAGRAM  MANIPULATION  OPTIONS 

*  20  LIST  OPTIONS 

*  21  FORM  OLTF  =  GTF  *  HTF  (IN  CASCADE) 

*  22  FORM  CLTF  =  ( GAIN*GTF ) / ( 1  +  GAIN#GTF#HTF ) 

*  23  FORM  CLTF  -  GAIN*OLTF  /(IF  GAIN*OLTF) 

*  24  FORM  CLTF  =  GTF  +  HTF  (IN  PARALLEL) 

*  25  G'TF(S)  8  HTF(S)  FROM  CONTINUOUS  STATE-SPACE  MODEL 

*  26  GTF(Z)  8  HTF(Z)  FROM  DISCRETE  STATE-SPACE  MODEL 

*  27  WRITE  AD JOINT ( S I -AMAT )  TO  FILE  ANSWER 

*  28  FIND  HTF  FROM  CLTF  8  GTF  (FOR  CLTF=GTF*HTF/ ( 1 +GTF*HTF ) ) 

*  29  FIND  HTF  FROM  CLTF  8  GTF  (FOR  CLTF=  GTF  /( 1+GTF*HTF)  ) 

OPTION  >  25 

GTF < S )  AND  HTF(S)  CALCULATED  FROM  STATE-SPACE 
TYPE:  GTF  OR  HTF  FOR  RESULTS 


-35- 


FORWARD-LOOP  TRANSFER  FUNCTION 


r* 


GK=  <  GNK/GDK  >=  1.213 


GTF(S)  NUMERATOR 


I 

GNPOLY ( I ) 

GZERO(I) 

1 

( 

1.213  )S** 

8 

< 

. 8688E-1 1 ) 

+ 

J< 

6.880 

) 

2 

( 

84.99  >S** 

7 

( 

. 8670E-1 1 ) 

+ 

J( 

-6.880 

) 

3 

( 

5418.  >S** 

6 

< 

-19.81  ) 

+ 

J( 

1 . 554 

) 

4 

( 

.2494E+06)S** 

5 

< 

-19.81  ) 

+ 

J( 

-1.554 

) 

5 

< 

.6119E+07)S** 

4 

( 

-15.21  ) 

+ 

J< 

2.205 

) 

6 

< 

.7856E+08)S** 

3 

< 

-15.21  ) 

+ 

J< 

-2.205 

) 

7 

( 

. 5698E+09 ) S#* 

2 

( 

- . 4079E-12 ) 

+ 

J< 

50.81 

) 

8 

( 

.3168EH0)S** 

1 

( 

- • 4079E-12 ) 

+ 

J( 

-50.81 

) 

9 

< 

. 1383E+11 ) 

GNK= 

1 . 

213 

GTF(S)  DENOMINATOR 


I 

GDPOLY ( I ) 

6P0LE ( I ) 

1 

( 

1.000  )S**10 

(  - • 1916E-09 ) 

+ 

J< 

3.464 

) 

n 

( 

70.04  )S** 

9 

(  - . 1 91 6E-09 ) 

+ 

J( 

-3.464 

) 

3 

< 

6017.  )S** 

8 

<  -19.81  ) 

+ 

J< 

-1.554 

) 

4 

( 

.3142E*06)S** 

7 

(  -19.81  ) 

+ 

J( 

1.554 

) 

5 

( 

,9563E+07)S** 

6 

(  -15.21  ) 

+ 

J< 

-2.205 

> 

6 

( 

.2149E+09)S** 

5 

(  -15.21  ) 

+ 

J( 

2.205 

) 

7 

< 

.3703E+10)S** 

4 

<  • 1402E-1 1 ) 

+ 

J( 

-21.71 

) 

8 

( 

.3978E+11 )S** 

3 

(  . 1398E-1 1 ) 

+ 

J( 

21.71 

) 

9 

( 

.2056E+12)S** 

2 

(  -.2216E-12) 

+ 

J( 

60.80 

) 

10 

( 

.4470E-I  12)S** 

1 

C  - . 221 6E-12 ) 

+ 

J( 

-60.80 

) 

11 

( 

» 1951E+1 3 ) 

GDK= 

1 . 

000 

OPTION  >  HTF 


FEEDBACK-LOOP  TRANSFER  FUNCTION 


HK=  (  HNK/HDK  >=  14.31 


O 


HTF(S)  NUMERATOR 


I 

HNPOLY(I) 

HZEROCI) 

1 

<  17.36  >S** 

9 

(  . 81 1 IE 

-02) 

+ 

J< 

0. 

) 

2 

<  1276.  )S** 

8 

(  -7.386 

) 

+ 

J( 

5.583 

) 

3 

<  ♦  1169E  +  06JS** 

7 

(  -7.386 

) 

+ 

J( 

-5.583 

) 

4 

(  .5326Ei07)S** 

6 

(  -1.689 

) 

+ 

J< 

-12.35 

) 

5 

(  .  1601E  +  09)S** 

5 

(  -1.689 

) 

+ 

J< 

12.35 

> 

6 

(  ,2489E+10)S** 

4 

(  -22.63 

) 

+ 

J< 

-22.62 

) 

7 

(  . 3020ET 11) S** 

3 

(  -22.63 

) 

+ 

J( 

22.62 

) 

8 

(  ,2253E+12)S** 

2 

(  -5.053 

) 

+ 

J< 

-62.70 

) 

9 

(  .9350Ei 12)S*» 

1 

(  -5.053 

) 

+ 

J< 

62.70 

) 

10 

(  - . 7598E+ 1 0 ) 

HNK  = 

17. 

36 

HTF  <  S )  DENOMINATOR 


I 

1 

< 

HDF'OLY  ( I ) 

1.213  )S** 

8 

HPOLE ( I ) 

(  . 8688E-1 1 )  +  J< 

6.880 

) 

2 

< 

84.99 

)S** 

7 

(  • 8670E-1 1 ) 

+ 

J( 

-6.880 

) 

3 

( 

5418. 

)S** 

6 

(  -19.81  ) 

+ 

J< 

1.554 

) 

4 

< 

.2494E+06)S** 

5 

<  -19.81  ) 

+ 

J< 

-1.554 

) 

5 

( 

.6119E-I07>S** 

4 

(  -15.21  ) 

+ 

J< 

2.205 

) 

6 

( 

,7856E+08)S** 

3 

(  -15.21  ) 

+ 

J< 

-2.205 

) 

7 

< 

.5698E+09)S** 

n 

(  - . 4079E- 12 ) 

+ 

J( 

50.81 

) 

8 

< 

. 3168E+10 ) S## 

1 

<  -.4079E-12) 

+ 

J( 

-50.81 

) 

OPTION  >  If 


OPTION  > 
OLTF 


OPEN-LOOP  TRANSFER  FUNCTION 


OLK=  GAIN#  <  OLNK/OLDK  )= 
GAIN=  1,000 


17.36 


OLTF(S)  NUMERATOR 


OLNPOLY(I) 


OLZERO ( I > 


1 

( 

21.07  >S**17 

(  .811  IE-02) 

+ 

J< 

0. 

) 

2 

< 

3024.  )5**16 

(  -7.386  ) 

+ 

J( 

5.583 

) 

3 

( 

,3444E+06)S**1 

5 

(  -7.386  ) 

J( 

-5.583 

) 

4 

( 

.2764E+08)S**14 

<  « 8840E-1 1 ) 

+ 

J( 

-6.880 

) 

5 

( 

.1705E+10)S**13 

(  • 8839E-1 1 ) 

+ 

J( 

6.880 

» 

6 

< 

.8381E+11)S**1 

o 

(  -15.21  ) 

+ 

J( 

2.205 

) 

7 

( 

.3269E»13)S**U 

(  -15.21  ) 

+ 

J< 

-2.205 

) 

8 

C 

.9882E+1 4)S**10 

<  -1.689  ) 

+ 

J( 

-12.35 

> 

9 

< 

.2274E+161S** 

9 

(  -1.689  ) 

+ 

J( 

12.35 

) 

10 

( 

.4007E+17)S** 

8 

(  -19.81  ) 

+ 

J( 

-1.554 

) 

11 

( 

.5513E+18)S** 

7 

<  -19.81  ) 

+ 

J( 

1.554 

) 

12 

( 

.5984E+19)S** 

6 

(  -22,63  ) 

+ 

J( 

-22.62 

) 

13 

( 

» 5C73E+20 )  S#* 

5 

(  -22.63  ) 

+ 

J( 

22.62 

) 

14 

( 

.3319E+21 >S** 

4 

(  - , 2142E-1 1 ) 

+ 

J( 

-50.81 

> 

15 

( 

» 1664E+22) S## 

3 

<  -.2J04E-11) 

J( 

50.81 

) 

16 

( 

.6074E+22)S## 

2 

(  -5.053  ) 

+ 

J< 

-62.70 

) 

17 

( 

, 1290E+23 )  S#* 

1 

(  -5.053  ) 

+ 

J> 

62.70 

) 

18 

< 

- , 1051E+21 ) 

OLNK= 

21 

.07 

OLTF  <  S  > -  DENOMINATOR 


OLDPOLY ( I ) 


OLPOLE ( I ) 


1 

<  1.213  )S**18 

(  - . 1 916E-09  > 

+ 

J( 

3.464 

> 

2 

(  170.0  )S**17 

(  - , 1916E-09 ) 

+ 

J( 

-3.464 

) 

3 

(  , 1867E+05 ) S#*l 6 

(  -15.21  ) 

+ 

J< 

-2.205 

) 

4 

(  .  1521E+07)S)K*15 

(  -15.21  ) 

J( 

2.205 

) 

5 

(  .9449E+08)S**14 

(  . 8546E- 1 1 ) 

+ 

J( 

6.880 

) 

6 

<  .4783EH0)S**13 

(  . 8546E-1 1 ) 

+ 

J< 

-6.880 

) 

7 

(  . 1958E+12>S**12 

(  -19.81  > 

+ 

J( 

1.554 

) 

8 

(  ,6350E-H3)S**11 

(  -19.81  ) 

+ 

J< 

-1.554 

> 

9 

(  .1641E+15)S**10 

(  -19.81  ) 

+ 

J( 

1.553 

) 

10 

(  . 3422E+16 ) S*#  9 

(  -19.81  ) 

+ 

J< 

-1.553 

) 

11 

(  .5714E+17)S**  8 

<  -15.21  ) 

+ 

J( 

2.206 

) 

12 

(  .7453E+18)S**  7 

(  -15.21  ) 

+ 

J( 

-2.206 

) 

13 

(  . 7428E+19)S#*  6 

<  •  1656E-1 1 ) 

J< 

21.71 

) 

14 

(  .5674E+20>S**  5 

(  , 1 656E-1 1 ) 

+ 

J( 

-21,71 

) 

15 

(  .3414E+21 )S*#  4 

(  - . 901 3E-12 ) 

+ 

J( 

50.81 

) 

16 

(  . 1609E+22 ) S##  3 

(  -.1101E-11) 

+ 

J( 

-50.81 

) 

17 

(  . 5371E+22 ) S**  2 

<  . 5889E- 1 2 ) 

+ 

J( 

60.80 

) 

18 

(  . 1236E+23JS**  1 

(  . 5089E- 1 2 ) 

+ 

J( 

-£>0 . 80 

) 

19 

(  . 2697E+23 ) 

OLDK= 

1 

.213 

r 


FACTORED  INPUT  OF  OLTF 

ENTER  NUM  &  DENOH  DEGREES  (OR  SOURCE):  9»10 
ENTER  NUMERATOR  CONSTANT:  21,07 
ENTER  EACH  ROOT— RE  r  IN 


\ 

:  'K 

i 

OLZERO( 
OLZERO ( 
OLZERO ( 

1)  = 
2 )  = 
3)  = 

.0031111 ,0 
-7. 386 >5. 583 
(  -7.386 

>  + 

J< 

-5.583 

)  ASSUMED 

OLZERO ( 
OLZERO ( 

4 )  = 
5  >  = 

-1 .688>  12.-35 
(  -1.689 

)  + 

J( 

-12.35 

)  ASSUMED 

' 

OLZERO ( 
OLZERO ( 

6 )  = 
7  )  = 

-22. 63 » 22. 62 
(  -22.63 

)  + 

J( 

-22.62 

)  ASSUMED 

j 

OLZERO ( 
OLZERO ( 

8)  = 
?)  = 

-5. 053 » 62. 7 
(  -5.053 

)  + 

J( 

-62.70 

)  ASSUMED 

OLTF  NUMERATOR  ( OLNPOLY ) 


OLTF  ZEROS  (OLZERO) 


r  i 

r-  '  l 

( 

21.07 

)S## 

9 

( 

.  81 1  IE 

-02) 

+ 

J( 

0. 

< 

1549. 

)S## 

a 

( 

-7.386 

) 

+ 

J( 

5.583 

*  ' 

( 

.1419E+06)S** 

7 

( 

-7.386 

) 

+ 

J( 

-5.583 

( 

.6464E+07)S** 

6 

( 

-1.689 

) 

+ 

J( 

12.35 

. 

( 

. 1944E+09)S** 

5 

( 

-1.689 

) 

+ 

J( 

-12.35 

> 

( 

.3021E+10)S** 

4 

( 

-22.63 

) 

+ 

J( 

22.62 

< 

• 3666E+1 1 ) S#* 

3 

( 

-22.63 

) 

+ 

J( 

-22.62 

< 

• 2735E+12 ) S## 

p 

( 

-5.053 

) 

+ 

J( 

62.70 

< 

. 1135E+13)S#* 

1 

( 

-5.053 

) 

+ 

J( 

-62.70 

<  - . 9221E+10 ) 


POLYNOMIAL  CQNSTANT= 


21.07 


ENTER  DENOMINATOR  CONSTANT:  1.IJ3 


ENTER  EACH  ROOT — RE  * IM 
OLPOLE  <  1>  =  0 » 3 . 464 


OLPOLE ( 

2 )  = 

(  0, 

)  + 

J( 

-3.464 

)  ASSUMED 

OLPOLE ( 
OLPOLE ( 

3 )  = 

4 )  = 

-f 9 i81 . 1 . 553 
(  -19.81 

>  + 

J( 

-1.553 

)  ASSUMED 

OLPOLE ( 
OLPOLE ( 

5 )  = 

6 )  = 

-15.21 »2. 206 
(  -15.21 

)  F 

J( 

-2.206 

)  ASSUMED 

OLPOLE ( 
OLPOLE ( 

7) ~ 

8)  = 

0»21 .71 
(  0. 

)  + 

J( 

-21.71 

)  ASSUMED 

OLPOLE ( 
OLPOLE ( 

9 )  = 
10)  = 

0 » 60 . 8 
(  0. 

)  + 

J( 

-60.80 

)  ASSUMED 

OLTF  DENOMINATOR  (OLDPOLY) 


OLTF  POLES  (OLPOLE) 


( 

1.213  )S**10 

( 

0. 

)  + 

J( 

3.464 

) 

( 

84.96  )S** 

9 

( 

0. 

)  + 

J( 

-3.464 

) 

( 

7298.  )S** 

8 

( 

-19.81 

)  + 

J( 

1.553 

) 

( 

.3810E+06)S** 

7 

( 

-19.81 

)  + 

J( 

-1.553 

) 

< 

.11 60E+08 ) S## 

6 

( 

-15.21 

)  + 

J( 

2.206 

) 

< 

« 2606E+09 ) S#* 

5 

( 

-15.21 

)  + 

J( 

-2.206 

> 

< 

.4491E+10)S## 

4 

( 

0, 

)  + 

J( 

21.71 

) 

( 

•  4824E+1 1 ) S** 

3 

( 

0. 

)  + 

J( 

-21.71 

) 

( 

.2493E+12)S## 

2 

( 

0. 

)  + 

J( 

60.80 

) 

( 

,5419E+12)S## 

1 

( 

0. 

>  + 

J( 

-60.80 

) 

( 

.2365E+13) 

POLYNOMIAL  CONSTANT= 

1.213 

GAIN=  1.0  OLK=  GAIN# ( OLNN/OLDK ) =  17.37015663644 

OPTION  >  AA  =  10  THESE  COMMANDS  SET  THE  PLOT  BOUNDARIES 

AA-  10.00000000 

-38- 


OPTION  >  BB=70 


UL<= 


70.00000000 


OPTION  >  DD--70 

DD~  -70.00 000000 

OPTION  >  PLOT t ON 

OPTION  >  ANSUERfON 

ANSWER  ON — OUTPUT  WILL  GO  TO  FILE — ANSWER 
OPTION  >  42 

ENTER  GAIN  OF  INTEREST  < GAIN ) f TOLERANCE  <GTOL)J>  If. 01 


TITLE  OF  THIS  PARTICULAR  PLOT — PRINTED  INSIDE  BOXJ 


— 

- ENTER 

TITLE 

(50  CHARACTERS 

MAX) - 

- 3 

ROOT 

LOCUS f  3M 

BEAM 

WT  =  DIAG<500 

500 

500 

2000) 

ROOT 

LOCUSf  3M 

BEAM 

WT  =  DI AG  <  500 

500 

500 

2000) 

OPTION  >  STOP 


LOCAL  FILE  — PLOT  — CONTAINS  CALCOMP  PLOT(S) 

LOCAL  FILE— ANSWER— CONTAINS  OUTPUT 

ALL  INFO  IN  TOTAL  HAS  BEEN  SAUED  IN  LOCAL  FILE-MEMORY. 
STOP 

9.632  CP  SECONDS  EXECUTION  TIME 
COMMAND-  ROUTE r PLOT f ST=CSB f DC=PT f  TI D*BB 

COMMAND-  ROUTE f ANSWER fDC=PRfST=CSBfTID=BBfFID=F 11 

COMMAND-  LOGOUT 


r > 


-39- 


o°pn--o-i°  n.rn  poo-  locus  using  option  4p 


:  00^ 

O-LOC'i7*  pt, 
1 -POL  E 


2-7  EPO 
7 -PEE AK  pt. 


4-  IH  AG  1*11  ?  V  AXIS 
li-SENSITWI  T  Y  &T. 


10  POLES 
y  = 


a  r 

o . 

0  . 

-19.110 
-1  G.M  0 
~  1  5.210 
-1  i  ."*10 


3.-40 
"3.4:  40 
1  .Or  70 
-1.  3  0 

2 .  7  r  s  n 
-2.2i.-b0 
PI  .710 
-2i. ;  10 
60. EDO 
-60.-00 


PECOS  \T 

y  =  .  3  in i ! 

y  =  -t .  33  o 

X  =  -  7 . 3  3  b  0 

X  =  -1.FS39 

X  =  - 1 . C  03  n 

X  =  -?P.670 

X  =  -IIo30 
X  =  -5.0-X0 

x  =  -r . 0"  3n 


0. 

G .  v  ?  3  0 
-D."ST0 

1?.  -so 

-i  c. o 
22. r  20 
-2  2.  ECO 
f>?.  ’0  0 
-b  2.730 


gain  CONSTANT  (olnx/qlok? = 


1  7  .3701  = 


GAIN  Oc  IN7ERrM  (GAIN)  = 


1.  ono 


SENSITIVITY  0P  I N  T  r  R  F  S  T  { QL  K)  = 


17.  V 


♦  OR-  GTOL  =  *  1  0  0  0 E -! 

♦  DR-  KT0.  =  .1737 


ROOTS  O-  T NT :  ? c >r 
y  =  -3.1’jn 

X  -  -3.13G0 

x  =  -  .  7  p  1  1 

y  =  -ll.  US 
X  *  -11  .116 

X  =  -  =>  .  ?  0  7 

X  -  -F.‘GT< 

y  -  i  .43  7 
X  r  .  1  t  A  -*  G 
X  =  .  1 1.  0  7  4 


-  ? .  S  L 1 5 

0. 

■1  3 «  r  7  •> 
1  3 .  •  7  ~ 
■?i.  <n  p 
?1 .407 

0. 

-S  c .  PRO 
6.2.24  = 


FUSION  OP 


:alcjlattom-f  "tL  j 

7  M  A  J  ! 


-9  c.  r 
-7  Q.o 


to  a  a t 

TO  3P  = 


">  G  A  N  C  H  TLUlMR  1 


■U  C’JLAT  I 


T  r  °  '  T  7  <r  =  .  r  0  G  0 


0. 

7 . 't  -  ■ »  - ;  •.  g 

:  ,  ‘  >  -4  0  03  0 

0  . 

.  0  '  fl  4  4 

-.3.0/7;- 

•• . '  3  3  r  ;> 

.  ■’9  o  ?r  i 

.  /  -  ? 1 1 

-3.2323  1  '*  ‘r 

*>.  *  -  ? 

147:;677 

1 .  0  l‘9b 

.  7 o  r  ?  j 

-5-.',M??7 

4 .  ’  •*.  V  ’  -  1 

7 

2  .  2?  06  9 

.  7  •  > .  i  7 

-7.7350010 

i  fl 

9.  ?«;■}  5654 

0  . 

. 79774 

’c-AUGM 

Mwr^  ? 

CU9tlL4Tr 

n»i  ;Tn  ~  T  ?  f_  = 

.5009 

3  ■<  1 911:19 

;  T  ro  <;!•"'  i 

4.000 

nm  mi 

_ 0  7 1  <  7  T430 

0TGT  ro  0=71319 

G?IM 

7  ETC 

0. 

-  7 . 4  r>  i  n  i  o  o 

3.-. s-t  o  ro  j 

0  . 

.  0  0  C  0  0 

-3.077V-?7 

-  2 . 5  1  ?  ?  “  '  a 

*>.  (m03?4 

.'•-9  0  2:1 

.70713 

>3.  U5-i 

"2 

4 « 1  4  3  5  b?  7 

1 .01445 

.70020 

-3.  mi?7? 

*  '  ♦  '  1  ^  1 

7.  ??'»?&!  f 

? .270:0 

•  7 1) .  7  0 

-7.33-- a  010 

-T/'nina 

9.25  3  666-4 

0  . 

.7977 + 

9PANCM 

rU.H9cn  3 

3U.CUL1TI  04 

i  >  T  E°  "T7E  = 

.800  3 

piruiir,  s t 

ro  q7r  _ 

4.000 

lo:j3  mi 

-.9-'')?  I  rt,09 

D1ST  TO  0P.I3I» 

GUM 

’5T  7 

-  19.31-3  000 

l .  *  p  3  n  o  c  c 

19.370730 

0  . 

.9?r 0* 

-19.-67  3‘  » 7 

0. 

19.473547 

.1299556-0*’ 

L  .  OCC  0  0 

-73.4?  3  ;  4  7 

0 . 

23.  47  35*7 

.  1744?  2  5-01 

1 . 0  0  r  0  0 

-  27  .  '4  7  ,3  4  7 

0. 

27. 4735*7 

.  c  6  7  2  5  4  E-  0 1 

1 . 0  c  m 

-31.473'!  47 

0. 

31.473  It? 

. 269939 

1.  ccc-oo 

-3?.-,?  39V 7 

0  . 

3V.47  3--7 

.“ 21912 

1.00m 

-7<».47  ?c47 

0. 

39. 4?  3  E-t  7 

.  °267 12 

1 .  0  0  r  0  D 

-41.2? 3 v*7 

0  . 

i  1.3 7  35  V 7 

1 .  0  2  *‘7  3 

i .  ocr  00 

-•4r  .  ?  7  ?  '?  4  7 

0  . 

45.233  347 

1.30914 

1 . 0  0  r  ■:  0 

-49. 27  4? 

0  . 

49.273547 

1 .05059 

1.0000° 

-5-3. 2?  3  447 

f» 

••  * 

53.7?3547 

1 .95770 

1.0cm 

-57.2??r4? 

0. 

57.2  33  547 

2.316^6 

t ,  0  c  0  0  0 

-  b  1 . 2  3  3 1‘  4  7 

0. 

61.2335 47 

2.63049 

1 .  c  c  0  0  0 

-65.27314  7 

0. 

05.273547 

2.94  777 

1  .  00O-) !) 

-69.23  3-4? 

n 

■  • 

69.2  33  54  7 

?.249:9 

1.00000 

-7  3.  ?3  3  '  4  *7 

0  . 

73.273597 

3.3*445  9 

1 . 0  u  c  j : 

-77.2? 3' 4? 

0. 

77 .733547 

3 . 3  3  2°  1 

l.OOC'G 

-81.23  31=47 

0. 

P1.23?5<*7 

4.11517 

1  .  J  0  r  0  n 

-30- . .? 7  3 •-  *  7 

n  . 

to.  273  54  7 

4  .  39  24  0 

1.  com 

-39.273943 

0  . 

69.273  6-4? 

*•  .66911 

1.00000 

-90.033!  ►? 

o. 

90. 03 ?5h7 

4,7 1916 

1.00m 

90'j-'in»  pv 

iJ  O 


•  <  1  -  ■>  '. 
..uuuir-1' 


•  3  1 1  1  1 1]  0  0  -  -  n  ?  n  . 


1.00(00 


a r t ;Tr'>  ■ 

MINTING  ST^  ^T7: 


anch 

7  r  =  .  l  c  c  o 

=  t*  .033 


LOC’J  5  ?-M. 

■  l’j.’nn11 
■lb.  Si  -*7  3  < 

h  U 

■  n .  i  7 1 7  5  4 

•  1  o  .  rr  »  *  * 6 

■  ?  .  5  i  3  '1  .  7 

■  7.P,;,  ,051 

■  1  . '  3  9  1  r  0  3 


.  103  IMAG 
7.  ?  9  5  .V*  C  0 

4.  '  ■  17  ^  ->  X 
°  r  r 

li  .f  4  =  7iv 
1  ?  .  jvtpo 

1*  .  •’I-Vi. 

1  ?  .  ’  T'  '  <  -  b 

i 7 . 7 : 3 '  o  c 


HIST  TO  0.1  31  N 
U  .75  9  l'»7 
1  7  .  ?  f  7:  r-  , 
17.-6  ♦  ■  ?-* 

IV.  -44  7 

1 7 .  -  1  7  t-  7  2 

lv.  Q0l*t-- 

17.77 
17.  4b496Q 


GAIN 

0  . 

.<66? 75 6-01 
.  3 1  C  3na 
,r3Si^; 

1  .  C  f  7  ?  6 

1 .  ? »  1  y  ? 

*•  .00796 

0  . 


7  FT  A 
. 96635 
. 9' <  ’a 
. oV673 
.  7  *-7  7  1 
.61!  1  " 
.  «  6  ■  5  3 
. 2°C  4  3 
.17t;3 


’ano 


3  UGULATf  on  ;•"-">  7T7p 
MINTING  GTE  =  SI7r 


L03JT 

■ls.mr  30 

■  1 6 . 5  a  7  ’  3  '3 
•  1 5 . 7  7  i  f  u 
■17.17  3  7  ? 

•10,7/?F7t; 

■7.733  -J  ■■  7 

■  7 .  3 "  ■  -» c  7  7 


LOO?  T  *<£G 
-7.  7  1.3  7  ‘  OC 
-  4 .  •  1  -  r’  7  7 

-0. T?1G- 7b 
-11  .  c  u  7  O  7 
-17  . -7 ' " 3  0 

-1  ♦  .  7m:o 

-1  ?.  7  1U-/.C 


.  b  0  3  0 
4.000 

01  Si  TO  ORIGIN 
13.  359  1  42 
17.7?7  964 

17.  A  "  4.  -»?4 

1 7 .  b  *  *■  9  6  7 
17. SI Zh?2 
Id.  0  9  4  3  ‘♦‘+ 
13.27  066 2 


GAIN 

0  . 

.  166726  6-01 
. ! 5C°! 6 
, r  94 123 
1 . 05778 
1 . 31- 1:  :j 
b. 0  0  7*5 


7  ETA 
. 935  5 : 
.961-79 

. 37t  75 
.74731 
.  611-1? 
.  46C-33 
.79040 


V-  ■ 

!■ 

i 

“1*6333030 

-1?  .  310  0  6  0 

1 2 .  h  6  4  9  5  0 

0  . 

.  17?  53 

’ 

i  1 

1 

np  anch 

OALC'JLATIOM  ;7-o  ~j7r  _ 

.eouo 

1 

| 

ANTING 

>Tr3  ^T7-  = 

4 . 0  3  D 

m 

LOOJT  9^1 

L  Of"  1  -  Tvtf.G 

GIST  TP  0  31 ‘-IN 

GAIN 

7  ETA 

-4 

0. 

?l .  7 1  m  0  c 

21.713000 

0  . 

. 00L00 

i 

-  3  .  G  '  I  .  6  7  1 

7t  .  4  .7-.1  7 

2  2.  207  -53 

.r  43106 

. 17? 05 

j 

-7.1 7'  3  7,7 

7t  .  731  ' 

27.  4 > 240 

1.06474 

.  317  4  9 

-10.  *307/1 

-'l  ,  r  -»  7  T  r  -> 

?7. 032  69  9 

1 .4  61'-.  0 

.  4  3  r  ?  9 

-14.  J  4  M37 

2  3  .  7  7  •:  1 7  u 

27. 1  57  206 

?. 22723 

.511  7” 

-17.70  3 1  ->  S 

7  +  .  7'7  '  7R 

30. 1?  r  ir»  7 

4.1^249 

.50790 

-71  .’7  747  1 

31. 16 0-61 

17. 5657 

.63153 

-  2  2 .  •  1 3  0  0  0  0 

7  ? . 6 70100 

3 1.9? 66 3 7 

0  . 

.70775 

OILO'JLUIC"  I7"-1  —  .oOOO 

^'lINTI*'-,  ~  T  r  7  91  7"  =  4 . 0 c  0 

LOO  J  5  =>cAl  .I^T  t  b  A  g  pis:  to  oeigtn 
-21.7  1*’ ■>0  0  21./  1  3  COO 


GAIN 

Q  . 


7  CTA 
.  0  0  r  0  0 

.1750& 


I 


I 


17  i  ) 1 1 !  j  i 

-  ->  j  .  '  -  t  i  ^  r. 

■  r .  i  s  o  7 

<■  .  1 

.  9  9  7 

j  j 

2  1 .  '  ’  ’ '  ’’I 

">]  7(. 

'  i  .  7 o  '-M 

1  ~  .  :>  ?  b  2 

.  t>  i  1 

3  7 

2  0  ]  o 

.  -MOO  c  c 

7 1 .  ^  3  f.  :•  i  3 

P  . 

.707 

21 

■'PAUO 

1  9 

*"*  A 

L  D 1  i  L  ' 

riot:  ' 

T"">  pI7f  _ 

.  P  C'  0  0 

INTI  N 

>T'  = 

«;t  = 

'•.000 

LOCJ~ 

TPM 

nnj- 

L'isr  rc  o-i i” in 

GAIN 

7 

FTP 

0  . 

'  ] 

. 6  10  ■>  n  o 

f.f’  .700  on  o 

P  . 

• 

0  ?C  G  0 

.12297 

L  ?  1 

4? 

.  7  7  ’  •  «•  0 

7  7  7  f  0  1 

1 . 09o39 

• 

0  01-07 

0. 

•'•>? 

.  ”  v»  *•'  t 

l-  2.  7  7UOE7 

l.upu7u 

• 

o  o  r ;  o 

7. 2*7  ^ 

'  "» 

r:  '-4 

C  4 . 5  3  f>  4  7  0 

7  .?'  93  3 

• 

o tr  r 7  2 

■5.0 

’  3  0 

D? 

,  -  7  n  -  0  0 

f-  2.  3  0  3  221 

0  . 

• 

op.p  ;<7 

”*P  A  MCH 

NU  *7  PE  3  10 

D  ALGOLS  T I  ON  3  T  r">  ’T?E  = 

.0000 

MINTING 

PJP5  r~y  7  ~  — 

7.  .00  0 

LOG  1 3  rAL 

_  a  ?  1 7  :  n  a  r. 

DIET  T  3  DPI  7 I M 

PAIN 

7  ETA- 

0. 

-07 , 1 1 p  t  n  o 

»  0.  ?  0  0  00  0 

C  . 

. OOCOO 

.122770  i  1 

92. 3’7301 

1 . 09o39 

.00197 

G. 

-O’  .771.7-8 

62.7  Ki»93 

1  .4,2474 

. cooo" 

3. 27  77  03? 

.  7  O'-  75  O 

64.  3Jtr<.7 0 

7 .24999 

.  CM  2? 

5. 057  3  00  0 

-3?  .70n"OG 

6  2.307231 

0  . 

.  C?  C  7  ? 

ciijpjj 


////  END  Of  LI5T  //// 


ROOT  LOCUS  FOR  3m  BEAM 


WEIGHTING 


_ _ . _  f 

-21.43  -14.29  -7.14 


ROOT  LOCUS  FOR  3m  BEAM 
WEIGHTING  =  DIAG  (500  500  500  500) 


-25.71 


-17.14 


-8.67 


Fi'furo  Hi.  !t»t  kxjus  diaqrams  har  hi  sc  wciqhtinqs 


J' 


-25.71  -17.14  -8.57 


ROOT  LOCUS  FOR  3m  BEAM 
WEIGHTING  =  01 AG (100  111) 


-47- 


-48- 


Figure  B5.  toot  locus  diagrams  for  fourth  state  weighting  variatio 
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ROOT  LOCUS  FOR  3m  BEAM 


Figure  B6.  Root  locus  diagrams  for  first  state  weighting  variation 
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Figure  B7.  Root  locus  diagrams  for  second  state  weighting  variation 


Figure  B8.  Root  locus  diagrams  for  tliird  state  weightina  variation 


Figure  B9.  toot  locus  diagrams  for  fourth  state  weighting  variation. 
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Figure  BIO.  Root  locus  diagrams  for  large  first  two  state  weightings 
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Fiqure  Bl2.  Time  response  for  weighting  =  DIAG(500  500  500  500) . 
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Figure  B14.  Tims  response  for  weighting  =  DIAG(100,000  500  500  500) 
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Figure  B18.  Time  response  for  weighting  =  DLAG(500  10,000  500  500) 
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Figure  BL9.  Time  response  for  weighting  =  DIAG(500  10,000  500  500) 
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E’igure  B21 .  Time  response  for  weighting  =  DIAG(500  500  2000  500) 


Figure  B23.  Time  response  for  weighting  =  DXAG{500  500  500  100). 
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Figure  B25.  Time  response  for  weighting  -  DIAG(500  500  500  100) 
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Fig  b26 .  Sensor  displacement  for  one  ze 


FOR  3M  BEAM  WITH  FEEDBACK  COMP:  (S+50) 
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B23 .  Sensor  displacement  for  one  zero  at  -50. 
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